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Abstract 

We study Dyson-Schwinger equations for propagators of Dirac fermions interacting with a mas- 
sive gauge boson in the ladder approximation. The equations have the form of the coupled nonlinear 
integral Fredholm equations of the second kind in the spacelike domain. The solutions in the time- 
like domain are completely defined by evaluations of integrals of the spacelike domain solutions. 
We solve the equations and analyze the behavior of solutions on the mass of the gauge boson, the 
coupling constant, and the ultraviolet cutoff. We find that there are at least two solutions for the 
fixed gauge boson mass, coupling, and the ultraviolet cutoff, thus there are at least two fermion 
families. The zero-node solution represents the heaviest Dirac fermion state, while the one-node 
solution is the lighter one. The mass gap between the two families is of the order of magnitude 
observed in nature. 
PACS numbers: 11.15.-q; 11.15.Ex 



1 



I. INTRODUCTION AND MOTIVATION 



There is a common belief in particle physics that the Higgs mechanism resolves the 
problem of generating the masses of particles, although the masses of gauge bosons and 
fermions are fixed by couplings to Higgs scalars that are completely free parameters. The 
observed pattern of fermion and gauge boson masses should stimulate a search for a better 
symmetry breaking mechanism. The example of such a symmetry breaking mechanism is 
the proposed principle of noncontractible space [lj]. 

It is assumed that the origin of the broken conformal, discrete, and gauge symmetries in 
particle physics is hidden in the character of the physical space, namely its noncontractibil- 
ity The appearance of the chirally asymmetric coupling of SU(2) gauge bosons to leptons 
and quarks, the appearance of the massive SU(2) gauge bosons and very heavy and very 
light Majorana neutrinos, the relation between fermion and gauge boson mixing angles, etc., 
are just consequences of mathematical consistency requirements of the proposed SU(3) con- 
formal unification scheme of strong and electroweak forces. The universal ultraviolet cutoff 
(minimal universal scale or distance) is fixed by the mass of weak gauge bosons 

2tt e 2 
A = M^, g w = e/ sin Q w , a e = — A = 321.3GeV. (1) 

V6g w 4?r 

It was shown that heavy Majorana neutrinos could be perfect candidate particles for 
cold dark matter because they are cosmologically stable owing to the absence of Higgs 
scalars in the theory. They are probably already indirectly observed in the center of our 
galaxy by atmospheric Cerenkov telescopes \&\. The universality of the minimal scale is 
confirmed by the Einstein-Cartan quantum theory of gravity that can resolve cosmological 
problems without inflaton scalar fields [5j]. The prediction of the Einstein-Cartan cosmology 
for the negative cosmological constant (or the negative contribution of torsion and zero 
cosmological constant) can explain the low power of the large-angle CMBR data by the 
integrated Sachs- Wolfe effect if the Hubble constant is small (or large). If the total angular 
momentum of the Universe is large at present then the Hubble constant could be also large 
explaining the anomalous large scale flows of the Universe [f|. The rotating Universe is 
a natural consequence of the Einstein-Cartan cosmology with spinninghot and cold dark 
matter particles {?]. This vorticity can be studied by CMBR (WMAP)[s| or by SDSS data 
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Having not only renormalizable theory 10|, but also ultraviolet finite gauge theory to 
describe the world of particle physics, one is faced with a possibility to resolve the com- 
plete structure of all Green functions of the theory studying the respective Dyson-Schwinger 
equations. We start this difficult task in the present paper by the study of Dyson-Schwinger 
equations for Dirac fermion propagators in the ladder approximation where fermions are 
coupled to one massive gauge boson. We give the respective equations and necessary algo- 
rithms to solve these equations in the next section, while the results and thorough analyses 
and discussion are given in the last section. 



II. EQUATIONS AND ALGORITHMS 

We assume that Dirac fermions couple chirally symmetric to one massive gauge boson by 
the standard form (M=gauge boson mass), thus we study the Abelian chirally symmetric 
version of the BY theory of the Ref. [fj: 



C = V D (ifi-g A)V D - + C g . f . 

+(d^* + igA^W® - igA^Q) - (Ym^d^d + h.c), (2) 
Cgj, = standard gauge fixing terms, ^> C D = charge conjugated $ £>, 
F^v = d^Ay - d v Ap, $ = v + ix, M — V2gv. 

The renormalizability and gauge invariance are ensured by the coupling of the Nambu- 
Goldstone boson to the gauge boson and the Majorana fermion as in We assume also 
that Nambu-Goldstone bosons carry lepton number as in the BY theory of Ref . [lj , thus only 
Majorana bare mass term is allowed. The symmetry breaking parameter v is in this model 
freeparameter and it is not fixed by the Wick's theorem as in the non-Abelian version BY 
of We study only equations of a Dirac fermion in this paper. 

The Dirac fermion propagator is defined as S' F (p) = [a(p) Jp — fiip)]^ 1 ■ It is advantageous 
;o write Dyson-Schwinger equations for fermion propagators in the ladder approximation 



in the Landau gauge. Then, the equations have the following form in the spacelike 
domain: 
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m = C f A2 dyK(x,y) J^l\ 2( y 
Jo yoL\y) + (J 2 {y) 

a( x ) = l-C [ dyL{x,y) — w^TTWi v 
Jo ya 2 (y) + (3 2 (y) 



(3) 



x 



' "~ • - <] K{x,y) = + ^~ 



p 2 , C = a g /n 
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3,. M 2 

2" 3s / x + y + M 2 + s'' 



L(x,y) 



yM 2 



s = [( x + y + M 2 ) 2 - Axy] 
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s(x + y + M 2 + s) 2 ' 

One has to add one more term in the timelike domain because of the existence of the 



branch point of kernels and a correct analytical continuation 10|, LL3j : 



«,) = c£\yK( X , v) ya ^ 2(y) 

+ 0(^5 - M)C J°^_ m ^AK {x ,«) - ? m . 



a(x) = l-C f dyL(x,y) 

Jo ya 2 {y) + (3 2 {y) 

- 0( v /z x -M)cJ dyAL(x,u) 



ya{y) 



ya 2 (y) + (3 2 (y)' 



(4) 



s = it, r = Axy -(x + y + M ) , AK = K(x, y, s) - K(x, y, s*), 



> AK(x,y) 
AL(x,y) = 



AL = L(x,y,s) - L(x,y,s*) 
3s M 2 (x + y + M 2 ) 



Axy 



Axys 



M 2 



Ayx 2 s 



[(x + y + M 2 ) 2 - 2xy] 



Now we shall describe in detail how we solve equations, while anayses of solutions are 
left for the final section. 

We solve equations in four steps. 
Step 1: 

The equations in the spacelike domain have the form of the coupled nonlinear Fredholm 
integral equations of the second kind [3] and we need the initial guess functions to proceed 
further. A good choice is a solution of the nonlinear equations for the vanishing gauge boson 
mass. In this case, the nonlinear integral equations are reduced to nonlinear differential 



equations 



This equation is solved by the Adams-Bashforth method. One can easily find initial 
conditions from the equation and using the rule of de l'Hopital at x = p 2 = 0: 



dB , . 3 C d 2 B . , 3C 2 

■(0)=-o^> -3^ 0) 



dx w 8 5(0)' dx 2V 7 32S 3 (0)' 

The asymptotic of the solution at large spacelike momenta gives us a condition to find 
the solution with an arbitrary number of nodes 

dB, a2n B(x = A 2 ) 
— (x = A 2 + ^ A2 l! = 0. 
dx v 11 A 2 

One can simply generate solutions from fixed initial conditions at the zero momentum 
and then search for momenta which fulfill above condition. As an example, for C = 0.7 and 
B(0) = lGeV, one obtains j-node solutions with the following Aj cutoffs: A = A.292GeV , 
Ax = 85.82Ge1/ and A 2 = 1715.87GW. 

Thus we can generate initial guess functions by (3 = B and a = 1. Assuming the 
existence of the fundamental cutoff A, one has to rescale all dimensional quantities of a 
certain dimension k by multiplication with (A/Aj) fc at the end of the calculation with the 
cutoff Aj. 

Step 2: 

Using the initial guess functions, we have to choose the method how to solve nonlinear 
integral equations. Between Nystrom, Galerkin or the collocation method 



121 . or Newton- 



like iterations [l3j, we decide to implement the collocation method. 

It is more comfortable to work with logarithmic variables in the spacelike domain, thus 
we change the variables to w — ln(l + x/B(0) 2 ). The (3 and a functions are approximated 
by Cebisev polynomials: 



1 N 

f(x) w --ci + c k T k -i{x) : -1 < x < +1. 
1 fc=i 
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Inserting (3 and a into the Fredholm equations, Eq.(3), we obtain a nonlinear algebraic 
system of equations for the set of coefficients \bj, a A in Eq. (6). This system of equations 
is solved by the modified Powell hybrid method 14]. 

We divide the integration region for a one-node solution into two segments < x < Aq 
and Aq < x < A 2 , and write the approximate solutions j3 and a with two separate Cebisev 
expansions, one for every segment. The accuracy of the approximation is improved, but the 
number of variables is twice compared with the zero-node case. Usually, we take 20 Cebisev 
polynomials in one approximation, thus 2 x 2 x 20 = 80 coefficients as variables for one- 
node solution, 40 for zero-node, 160 for two-node solution, etc. We define the discretized 
sequence of the corresponding logarithmic variables vj\ in each segment that is distributed 
homogeneously, except for a denser distribution of points at the joint of two segments. The 
corresponding number of nonlinear algebraic equations for coefficients hi, aj can now be 
formed: 

F M ee P( Wl )-CB(0) 2 rdw'e w 'K{x{w t \y{w')) J^gM Q 

Jo y(w')a 2 (w') + (3 2 (w') 

F 2ji ee a( Wl )-l + CB(0r rdw'e w 'L{x{ Wl \y{w')) "^"^1 = 0, 

Jo y(w')a 2 (w') + (3 2 (w') 

i = l,...,n, (6) 

n _ n 2w 
P( w ) = J2 hTk-i(w), a(w) = y^a k T k - X {w), w = 1, 

fc=l fc=l w a 

w A = ln(l + A 2 /B{0) 2 ), x{w) = B(0)\e w - 1), y{w) = B(0)\e w - 1). 

To this system, we apply the modified Powell hybrid method and verify the result of the 
computation. 
Step 3: 

If we obtain a solution even under small tolerance, one has to verify it with a very large 
number of arguments w (more than 1000) for both f3 and a and check the errors. One can 
expect larger errors for higher-node solutions. 

The final check whether our solution is a real solution or only some local minimum must 
be performed by solving the system once more, but now with a more precise approximation 
scheme, namely, the piecewise cubic spline method. 



The variables are now values of (3 and a functions evaluated for 200 arguments Wi and we 
form equations as was done previously for Cebisev coefficients. Thus, we solve an algebraic 
system with 400 variables using the modified Powell hybrid method. Note that for every 
solution ((3, a) there is also a solution (—(3, a). If the procedure does not diverge or does not 
end with a trivial solution, we have a very accurate solution in the spacelike domain that 
should be checked explicitly once more. 

Step 4: 

It is obvious that the relations in Eq.(4) are not equations for the timelike domain, but 
just formulas for evaluations of (3 and a in the timelike domain, knowing their solutions in 
the spacelike domain. 

Namely, let us define the following intervals: Iq = {x\0 < x < A 2 } and Ij = Iq{J{x\0 > 
x > — (j x M) 2 }. From the boundaries of integrals in Eq.(4) one can conclude that the 
timelike part of l\ is defined by integration over 7 , similarly, the timelike part of I 2 is 
defined by integration over I ± , etc. The procedure can be continued to arbitrary I k . 

Note that the integrals over timelike domains are very well defined by the endpoint 
singularity integration. The endpoint singularity at y — —{^/x — M) 2 of kernels AK and 
AL is integrable. 

If the propagator reaches the mass singularity for a certain x = —rrij G h, then at 
x = —(m,f + M) 2 G Ik+i the functions (3 and a diverge because the endpoint singularity 
becomes nonintegrable. However, the mass function /i(x) = (3{x)/a(x) is even at this point 
well defined as the limes of the quotient oo/oo. Anyhow, the numerics in the interval I^+i 
in the vicinity of x = — (m/ + M) 2 is becoming difficult, so we stop our calculation at Ik- 

We perform integrations using the Gauss-Kronrod and Clenshaw-Curtis methods which 
we have checked using the slow and accurate Monte Carlo method. The functions (3 and 
at are approximated in the timelike part of Ik by piecewise cubic spline with 800 values for 
each function. 

This completes our procedure of solving Dyson- Schwinger equations for fermion propa- 
gators. The next section is devoted to results and physical implications. 
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FIG. 1: This figure depicts the mass-function /j, (solid line) and the j3- function (crosses) in the 
spacclike domain for a zero-node solution and parameters C = 0.7, M = 2>.f>GeV, A = 321.3Gey. 
"w" is a logarithmic variable defined by B(0) = lGeV. 

III. RESULTS AND DISCUSSION 

The results presented in this paper are only those that have passed our procedure in four 
steps. The fact that the procedure fails in certain cases does not mean that the equations 
do not necessarily have solutions. The proof of the existence of any kind of solutions is not 
attempted in this paper. 

Our experience with these equations tells us that our procedure fails to find solutions 
for higher gauge boson mass, one- node solutions for coupling close to the critical C — |, 
and two- or higher-node solutions for any coupling. However, we find enough solutions to 
make relevant physical conclusions and suggestions for the improvement of algorithms. The 
attempt to find a solution by perturbing parameters, such as coupling or boson mass, is 
usually not successful. 

The reader can visualize solutions with zero and one node in Figs. 1-4, where the mass 
function is defined as \i = (3 /a and rrif = fermion mass. 

That zero-node solutions are heavier than one-node solutions can be seen from Figs. 1-5. 
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FIG. 2: This figure depicts the mass- function fi in the timelike domain for a zero- node solution 
and parameters C = 0.7, M = 3.5GeF, A = 321.3GeF. 
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FIG. 3: The mass- function fi (solid line) and the /3-function (crosses) in the spacelike domain for 
a one-node solution and parameters C = 0.7, M = 3.5GeV, A = 321.3Gey. "w" is a logarithmic 
variable defined by B(0) = lGeV. 
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FIG. 4: The mass-function n(x) (solid line) and sf^x (crosses) plotted in the timelike domain for 
a one-node solution and parameters C = 0.7, M = 3.5GeV, A = 321.3GeV. 
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FIG. 5: The quotient of the mass- functions fj, of zero-node and one-node solutions in the timelike 
domain for parameters C = 0.75, M = 3.0GeV, A = 32l.3GeV . 
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FIG. 6: The quotient of the mass-functions jj, of zero-node solutions for C = 0.7 and C = 0.5 in 
the timelike domain and parameters M = 4.5GeV, A = 321.3GeV . 



M(GeV) 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


m f (GeV) 




4.932 


2.528 


1.912 


1.310 


0.642 



TABLE I: The dependence of m/ on the gauge boson mass M for a one-node solution with 
A = 321.3GeF and C = 0.7; sign "-" denotes the absence of the mass singularity 

Equations (3) have solutions only in the strong coupling regime C = a g /n > |. Verifying 
the asymptotics of the zero-node solutions in the timelike domain, one can conclude on 
the absence of the mass singularity for couplings or gauge-boson masses when we have 
nontrivial solutions. This is a consequence of the strong coupling between fermions and the 
gauge boson. Usually, the one- node solutions have mass singularity (see Tables 1-3), unless 
one reaches a very strong coupling regime. 



c 


0.65 


0.675 


0.70 


0.725 


0.75 


m f {GeV) 


0.351 


1.401 


2.528 


5.833 





TABLE II: The dependence of rrif on the coupling constant C for a one-node solution with 
A = 321.3GeF and M = Z.QGeV 
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A(GeV) 


214.2 


241.0 


275.4 


321.3 


482.0 


m f {GeV) 


0.428 


0.983 


1.639 


2.528 





TABLE III: The behavior of mj on the fundamental cutoff A for a one-node solution with M 
3GeV and C = 0.7 



M(GeV) 


6.0 


8.0 


10.0 


12.0 


14.0 


16.0 


18.0 


m f (GeV) 




22.43 


15.27 


23.90 


9.84 


7.24 


4.44 



TABLE IV: The behavior of mf on the gauge boson mass M for a zero-node solution with 
A = 321.3GeV and C = 0.5 in the "model" case a = 1 

The "model" case when a = 1 is studied for zero- node solutions with results in Table 4, 
where one can see the appearence of the mass singularity and its behavior on the gauge boson 
mass. The same "model" is studied in [l5] but it is deficient because of various reasons: (1) 
it is not gauge invariant within a ladder approximation, (2) nontrivial solutions emerge only 
in the strong coupling regime, so a function affects crucially the mass function and cannot be 
set to 1, (3) according to the discussion in Step 4 of the description of the solving-procedure 
it is evident that the mass function contains spurious unphysical singularity in the segment 
Jfc + i in the timelike domain if we put a = 1. 

The behavior of the mass of the one-node solutions on the gauge boson mass can be read 
in Table 1. Fermion masses are very sensitive to the boson mass and they are larger for 
smaller boson masses. This is the expected feature, because the fermion mass singularity 



disappears for the vanishing boson mass 



Larger coupling implies larger fermion self-energy, thus also larger fermion masses, as one 
can read from Table 2. By rescaling of cutoffs and masses in Table 1, in Table 3 we show 
the dependence of fermion masses on the fundamental cutoff for the fixed boson mass. The 
guess that larger cutoff means larger self-energy and consequently larger masses is completely 



c 


0.65 


0.65 


0.70 


0.70 


0.75 


0.75 


n 





1 





1 





1 


buv 


-1.003 


-1.0015 


-1.0021 


-1.00937 


-1.0014 


-1.0076 



TABLE V: The UV asymptotic index bjjv f° r various couplings C, M=3 GeV and A = 321.3GeV; 
n=number of nodes. 
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confirmed. 

Similar dependencies are valid for zero-node fermion states. We depict some comparisons 
of mass functions in Fig. 6. 

Let us define infrared (IR) and ultraviolet (UV) asymptotics of the mass function: 



. . ,.l , , xdu(x)/dx , , xdu(x)/dx 

li(p 2 ) ~ p 2 : b IR = hm —?y+ — , b uv = hm \ " — . 

A direct inspection of our solutions leads to the following conclusions: (1) bm = for 
M = and M ^ 0, (2) buv = ~ 1 f° r M = and there is a small deviation from -1 for 
M/0 (see Table 5). 

We can now make final conclusions about solutions of Dyson- Schwinger equations. We 
have shown that even the most simple bootstrap system studied in this paper has at least two 
solutions for the fixed coupling, the gauge boson mass, and the cutoff. The heaviest solution 
is the zero-node solution (belonging to the third fermion family in the standard particle 
physics classification), while the lighter one is the one- node solution (belonging to the second 
fermion family). The mass gap between these two solutions is of the order of magnitude 
observed in nature (Figs. 1-5). We need more sophisticated approximation schemes for 
functions to search for higher-node solutions, thus to find the first-family member fermion. 
The observed behavior of the mass functions on the parameters fulfills our expectations. 

One should consider our work as the starting attempt to solve the family replication 
problem which is not soluble with the Higgs mechanism. Two basic scenarios are possible: 
(1) there are only three physical solutions in the spacelike domain, (2) there are more than 
three solutions in the spacelike domain but only three physically acceptable also within the 
timelike domain. Future studies will answer which scenario will prevail. 

Before turning to a more complicated electroweak system, we are expecting to see the 
verification of the principle of noncontractibility by the LHC The nonresonant enhancement 
of the QCD amplitudes at the weak scale has been reported by the TeVatron [h]], especially 
measuring and estimating the quotient of cross sections at two center of mass energies 



171 ]. High-luminosity measurements at the LHC could determine the ultraviolet cutoff 
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to high accuracy. 
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